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MODEL STARS FOR THE MODELLING OF GALAXIES: 
q-ENHANCEMENT IN STELLAR POPULATIONS MODELS 

P. Coelho 1 
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ABSTRACT 

Stellar population (SP) models are an essential tool to understand the observations of galaxies and clusters. 
One of the main ingredients of a SP model is a library of stellar spectra, and both empirical and theoretical 
libraries can been used for this purpose. Here I will start by giving a short overview of the pros and cons of 
using theoretical libraries, i.e. model stars, to produce our galaxy models. Then I will address the question 
on how theoretical libraries can be used to model stellar populations, in particular to explore the effect of 
a-enhancement on spectral observablcs. 
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1. INTRODUCTION 

To understand how galaxies form and evolve is 
one of the long-standing questions in astronomy. 
One of the ways to address that question is to study 
the stellar content of the galaxies through the use of 
stellar population (SP) synthesis techniques, i.e., the 
modelling of the spectral energy distributio n (SEP) 
emitt ed by evolving stellar populations (see iTinslev 
1980). An evolutionary SP model has two main in- 
gredients: a set of stellar evolutionary models (tracks 
or isochrones) that will predict how the stars are dis- 
tributed in the HR diagram, and a library of stellar 
observables (e.g. colours, spectra, spectral indices, 
etc.) that coupled to the evolutionary models, will 
be used to predict the colours or spectra of a stellar 
population. By computing models for several choices 
of SP parameters (age t, metallicity Z, initial mass 
function IMF, etc.), it is possible to produce a grid 
of SP models that can be used to study a variety of 
systems, from early-type galaxies and spiral bulges 
to star forming galaxies at different redshifts. 

A keystone to extract the star formation history 
SFH of unresolved stellar systems is through the 
analysis of their chemical abundance pattern. The 
chemical pattern is a tracer of the history of star for- 
mation because different elements are produced at 
different timcscalcs during the evolution of a galaxy 
(e.g. Greggio & Rcnzini 1983, Matteucci & Grcggio 
1986). In fact, the a-elements over Fe abundances 
([a/Fe]) is largely used to constraint the formation 
timescale of a galaxy, since the a-elements (O, Ne, 
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Mg, Si, S, Ca, and Ti) are released early in the evo- 
lution of the galaxy by Type II supernovae (SNe II) , 
while Fe is mainly produced by Type la supernovae 
(SNe la) on longer time-scales. Galaxies that had 
undergone rapid star formation present [a/Fe] val- 
ues larger than the ratios found in the Milky Way 
disc, reason why SP models built with stars of the 
solar neighbourhood cannot reproduce the locus of 
more massive galaxies in Mg versus Fe indices plots 
(e.g. Worthev et al.lll992[ ). Therefore, large efforts 
have being made by different groups in computing 
SP models with different [a/Fe]. I will focus on the 
role that the theoretical stellar spectra computations 
have in modelling populations with different [a/Fe]. 

2. THEORETICAL STELLAR LIBRARIES 

The library of stellar spectra used in the SP mod- 
els can be either empirical or theoretical, a choice 
which is subject to debate. There are plent y of mod 



els in literature using both empirical (e.g. Vazdckis 



1999. the visible range of iBruzual fe Charlotl 12003 



and _ 

' e.g. UV and IR range of 



Le Borgne et al.l 12004) or theoretical libraries 



Bruzual fe Charlotl 120031 



Delgado et all 120051 iMarastonl 120051 . IZhang et al 



20051 ). An empirical library is based on observed 



stellar spectra covering as much as possible the pa- 
rameters T e ff, log g, Z. It is not a simple task to as- 
semble a library that simultaneously features high 
S/N, good flux calibration, large wavelength cover- 
age, high-resolution and accurately derived stellar 
para meters, but great improvem e nts have been made 
(e.g. iPrugniel fe Soubiranl |200lt IValdes et all 12004 ; 
Sanchez-Blazquez et al.l 20061) . The major drawback 



of an empirical library is that such high quality ob- 
servations are limited to the closest stars, and there- 
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fore the coverage of the HR diagram and the stellar 
abundances are biased towards the solar neighbour- 
hood. 

A theoretical stellar library (also called synthetic 
stellar library) is based on model atmospheres and 
atomic and molecular line list. A model atmosphere 
is the run of temperature, gas, electron and radia- 
tion pressure, convective velocity and flux, and more 
generally, of all relevant quantities as a function of 
some depth variable (geometrical, or optical depth 
at some special frequency, or column mass). The 
synthetic spectrum, or flux distribution is the emer- 
gent flux computed based on a model atmosphere, 
and is required for comparison with observations. 
Theoretical libraries have the advantadge of covering 
the parameter space in T e ff, logg, and abundances 
at will. Moreover, a synthetic star has very-well 
defined atmospheric parameters and infinite S/N, 
and covers a larger wavelength coverage with (pos- 
sibly) higher resolution than a single observed spec- 
trum. To compute a large synthetic library can be 
demanding in terms of computational time, but it 
is usually feasible. There is a caveat though: be- 
ing based on our knowledge on stellar atmospheres 
and databases of atomic and molecular transitions, 
those libraries are limited by the approximations and 
(in)accuracies of their underlying models, and in fact 
we are not able to reproduce accurately all spec- 
tral types (see e.g. iGustafsson et al. 2007 ). There 
is a lot to be done yet in terms of including in a 
realistic way effects of 3D hydrodynamics, asphcric- 
ity, N-LTE, winds, non-radiative heating, chromo- 
spheric contribution, etc. Moreover, the databases of 
atomic and molecular transitions provide relatively 
few lines with highly accurate oscillator strengths 
and broadening parameters, besides being often in - 
complete (see the progress report by lKuruc3 2006 ). 
An additional point of confusion for the SP mod- 
eller (and user) is that theoretical libraries can be 
categorised into two groups: 

• libraries computed with a good treatment 
of the line-blanketing, therefore being good 
for spectrophotometri c predictions and low- 
resolution studies ( e.g. Castelli fc Kurucz 20031 



Gustafsson et al 
20051) : 



2003T 



Brott fc Hauschildt 



• libraries which arc computed with shorter, 
fine-tuned, empirically calibrated atomic and 
molecular line lists, being the ones more 
appropriate fo r high - resolution studies (e.g. 
Peterson et al.l l2001t Rodriguez-Merino et al 



are not suitable to predict colours due to miss- 
ing line-blanketing. 

The ability of some of the recent theoretical 
libraries in reproducing obser v ation s was recently 
assessed in iMartins fe Coelhol ( 2007 ). Concerning 
the libraries for low-resolution studies, broad-band 
colours predictions by three libraries were compared 
to an empirical UBVRIJHK calibration. Models 
can reproduce with reasonable accuracy the stellar 
colours for a fair interval in effective temperatures 
and gravities, but there are some problems with 
U-B and B-V colours and very cool stars in gen- 
eral (V-K > 3). The results for the B-V and V-I 
colours are illustrated in Fig. 1. As for the libraries 
aimed at high-resolution studies, IMartins fc Coelhd 
(|2007l ) analysed the performance of three libraries 
by comparing their predictions for spectral indices 
to measurements given by empirical libraries. Fig. 2 
shows the results for three indices against the empir- 



ical library MILES (|Sanchez-Blazquez et all l2006f) 



In general the libraries present similar behaviours 
and systematic deviations. In particular, the lists 
of atomic and molecular lines need further improve- 
ment, specially in the blue region of the spectrum 
and for the stars with T cS < 4500K. 

3. INCLUDING a-ENHANCEMENT EFFECTS 
IN STELLAR POPULATION MODELS 

As mentioned in the Introduction, the spectra 
of the stars carry the chemical signatures that are 
imprints of the SFH of the galaxy they belong to. 
Hence, and despite the current limitations of the the- 
oretical libraries as mentioned in § 2, we cannot rely 
solely on the empirical libraries if we intend to repro- 
duce the spectra of galaxies which have undergone a 
star formation different from that of the solar neigh- 
bourhood. In fact, there is more than one way of 
computing SP models with variable [a/Fe] values. 

3.1. The first SP models with non-solar patterns 
The first attempts to model SPs bey ond the 



solar-s caled pattern date back to mid-90s. iBarbuv 
(|1994l ) used synthetic stellar spectra to predict the 
Mg2 index for ages and abu ndances typi c al of those 
of bulge globular clusters. I Weiss et al.l (1995), us- 
ing observed spectra of bulge stars and for the first 
time including a-enhanced evolutionary tracks for 
a high mctallicity population, computed models for 
Mg2 an d <Fe> indices. Alt ernatively to the use of 
spectra, iBorges et al. ( 1995 ) adopted empirical fit- 
ting function^ with explicit dependence on [Mg/Fe] 



20051; ICoelho et al.ll2005l) . But in general thev 



2 The fitting functions describe how stellar spectral indices 
vary as a function of the stellar parameters T e g, logg, [Fe/H] 
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Fig. 1. Comparison between the colours predicted by 
synthetic fl ux distributions a nd th e colour-temperature 
relation by IWorthev fc Led (|2006l ). for stars represen- 
tative of a 10 Gyr population. Red-diamo nds cor- 
respond to ATLAS9 dCasteUi fc Kurucj 120031). green- 
triangles to PHOENIX (jBrott fc Hauschildd l2005h and 
the b lue squares to MARCS models (|Gustafsson et all 
2003). Filled and open symbols represent dwarf and gi- 
ant stars respectively. Filled black-circles are the values 
expected from the empirical relation. The bottom of each 
colour plot shows the residuals (difference between mod- 
els predictions and empirical calibration; the thin black 
vertical lines are the error bars of the empiric al calibra- 
tion). Adapted from lMartins fc Coelhol (|2007h . 

or [Na/Fe] to model SP indices with variable Fe, 
Mg and Na abundances. Other authors further 
explored the co nclusions and methods of the pio 



neer e ff orts, e.g., Idiart et al. ( 1996t h Vazdekis et al 



(|l997T l: iTantalo et all (|l998l ). Those first models 



were limited to a few indices or to a restricted range 
of SP parameters. 

3.2. Fitting functions & Sensitivity tables 

The study by iTripicco fc Belli (|l995T l opened a 
way to produce non-solar-scaled models for all the 
Lick/IDS indices. Based on synthetic spectra cal- 
culations, they quantified how each of the Lick/IDS 
indices changes with variations of individual chem- 
ical elements. These sensitivity tables (also called 
response functions), were computed for three com- 
binations of T e ff and log g, which correspond to a 
main-sequence dwarf, a turn-off star and a red-giant 
star of a 5 Gyr-old stellar population. 



or [Z/H] and sometimes, abundance [X/Fe] of another clement 
relative to iron. 
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Fig. 2. Comparison between th e predictions of thr ee 
synthetic libraries in the y-axis ([Martins et ai1l2005l in 
the lef t-hand column; iMunari et al 1 120051 in the middle 
colum; ICoelho et al.ll2005l in the right-hand column) and 
indices measured in the MILES empirical library in the 
a;-axis, for three spectral indices indicated in the figure. 
Blue squares are stars with T e ff< 7000K, green diamonds 
for 4500K < T eS < 7000K and red circles for T cS < 4500K 
(filled and open symbols are dwarf and giant stars re- 
spectively). The solid line is the one-to-one relation. 
The thick black symbols represent a Sun-like dwarf (dia- 
mond) and an Arc t urus-l ike giant (cross) . Adapted from 
iMartins fc Coelhol (|2007f l. 

B y combining e mpiri cal and theoretical informa- 
tion, Trager et al. (|200fj| ) proposed a method to use 
the sensitivity tables by Tripicco fc Bell ( 19951 ) to 
ap ply corrections t o the empirical fitting functions 



by IWorthev et alj (|1994l ). in order to produce SP 
model indices with arbitrary compositions H The 
method was extended by iThomas et al.l (|2003r l 
who published a large set of SP models for all 
Lick/IDS indices and variable values of [er/Fe], 
[a/Ca] and [a/N]. Until the present time, this 
technique is the most widely used to model SPs 
indices with variable abundances patterns, and 
both model indices and ingredients (sensitivity 
tables and fitting functio n s) are c onstantly updated 



Houdashelt et al. 2002; Tantalo fc Chiosi 



2004 iKorn et all l200l iLee fc Worthevl 1 2005 



Annibali et al. 2007: Martin-Hernandez et al 



2007 



Schiavonll2007l ) 

Important 
edge of the 



advancements in our knowl- 
chemical composition of galax- 



3 Alternatively, the oretical fitting func tions with explicit 
dependence on a/Fc l|Barbuv et alJ | 2003ft can be employed 
to co mpute o-enhanced models iMendes de Oliveira et all 
120051 1. but the theoretical fitting functions were limited to 
a few Lick/IDS indices. 
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ies have been achieve d through the use of 



those models (e.g. Kuntschner et al. 




2001; 


Proctor & Sansoml 2002: Proctor et al. 


2004; 


Thomas et al. 2005: Kelson et al. 20061 Smith et al. 


20061: de la Rosa et al. 20071). Nevertheless this 



approach is subject to some limitations: 

• the responses of three synthetic stars are used 
to correct all stars in an SP model and the im- 
pact of this approximation in the accuracy of 
the model predictions is uncertain; 

• only Lick/IDS indices can be modelled. Hence 
full-spect r um hi gh -resolution models s uch a s 

Bruzual fc Charlotl (|2003f k 



Vazdekisl (Il999h 



Le Borgne et al 



remain constrained to 
the solar neighbourhood chemical pattern, ham- 
pering the use of full-spectrum fitting tech- 
niques, colours or non-Lick indices in interpret- 
ing observations. 

• the effect of the abundance variations on the 
isochroncs is often not included in a consistent 
way. 

3.3. Fully theoretical approach 

Recent progress in the modelling of high- 
resolution stellar spectra opened the door to a 
new kind of models, in which the effects of abun- 
dance variations can be studied at any wavelength. 
This is enabled by the publication of several li- 
braries of theoretical, high-resolution stellar spectra 
for bo th solar-scaled and a-enhanced chemical mix- 
tures dBrott fc Hauschiidll2005l ; ICoelho et all 120051: 



Munari et al 



20051) . 

The first fully theoretical high-resolution SP 
m odels for an a-enha nced mixture were presented 
in Coelho et al.l ( 2007 ). who used an improved y er- 
sion of the stellar library in ICoelho et all (j2005h to 



compute SP models for solar-scaled and a-enhanced 
compositions for three values of iron abundance 
[Fe/H] and ages from 3 to 14 Gyrs. T hese models 



empl oy newly computed stellar tracks (jWeiss et al 
120071 ) with the same abundances as the stellar li- 
brary. For the a-enhanced mixture it is adopted a 
flat-enhancement of 0.4 dex for the abundance ratios 
of all the classical a elements (O, Mg, Si, S, Ca and 
Ti). The impact of the spectral and evolutionary 
effects is illustrated on the Figs. 3 and 4 for broad- 
band colours and spectral indices, respectively. The 
modelling of colours requires that both effects are 
taken into account, as they often have comparable 
effects. As for the spectral indices, the spectral ef- 
fect is overall the dominant one, but the evolutionary 




Fig. 3. Differences between the a-enhanced predictions 
and the solar-scaled ones, in units of the solar-scaled 
value, for broad-band colours shown on z-axis. Each 
row corresponds to a different [Fe/H] value, indicated 
in the figure. Left- and righ-hand columns show the pre- 
dictions for 4 Gyr and 12 Gyr SSPs respectively. The 
green-triangles and red-diamonds lines correspond to the 
evolutionary effect and spectral effects respectively. The 
black-circles line are the pred ictions when both e ffects 
are considered together. From lCoelho et all (|2007l ) . 

effect can be non-negligible (e.g. Balmer indices and 
near-lR indices). 

Another theoretic al study a t high- resolution is 



being developed by iLee et al" ( 2007 . Lee et al 



in prep.), based on th e stellar evolutionary tracks 
bv botter et all (|2007h and with the flexibility to 
explore the effect of abundance variations in an 
element-by-element basis. 

The fully-theoretical method has the advantage 
of providing a larger coverage in wavelength at higher 
resolution than the methods in § 3.2, and of pro- 
viding more accurate predictions given that the ef- 
fect of the a-enhancement on the stellar evolutionary 
tracks and spectral library is included in a consis- 
tent way. The caveat of this approach is that given 
the limitations of the synthetic libraries mentioned 
in § 2, these models are more affected by zero-point 
problems than semi-empirical methods. Hence they 
are not a straightforward replacement for the models 
based on empirical libraries, being more suitable to 
differential analysis 0. 



4 There are interesting work on deriving the stellar popula- 
tion parameters in a differential way, avoiding zer o-point prob- 
lems that affect both mode ls and data. See e.g. iNelan et al.l 
j2005l) ; lKelson et al.l 1120061 ). 
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Fig. 4. Differences between the a-enhanced predictions 
and the solar-scaled ones, in units of the solar-scaled 
value, for several spectral indices shown on z-axis. Each 
row corresponds to a different [Fe/H] indicated in the 
figure. Left- and right-hand columns show the predic- 
tions for 4 Gyr and 12 Gyr SSPs, respectively. The 
green-triangles and red-diamonds lines correspond to the 
evolutionary effect and spectral effects respectively. The 
black-circles line are the pred ictions when both e ffects 
are considered together. From lCoelho et all (|2007l ). 

3.4. Spectral corrections 

There is ongoing work that aim at combining 
the versatility of the method in § 3.3 in terms of 
exploring wider wavelength ranges at higher resolu- 
tion, and the accuracy of SP models based on em- 
pirical libraries. Similarly to using response func- 
tions to correct indices given by fitting functions, 
theoretical stellar libraries can be used to differen- 
tially correct SP models based on empirical libraries. 
To my knowled ge, there are two grou ps working on 



lo my knowled ge, there are two grou ps working on 
this approach. iPrugniel et al. (120071) recent ly used 
the theoretical library bv lCoelho et alj (|2005l ) to dif- 
ferentially correct the stars of the em pirical library 
ELODIE (|Prugniel fc Soubiranl [ioCJlh . Two 'semi- 
empirical' libraries were produced with [Mg/Fe] = 
0.0 and +0.4, and intermediate values were linearly 
interpolated as a function of the mass ratio of Mg to 
Fe. This semi-empirical library can then be used 
to build SP models with variable [Mg/Fe]. Al- 
ternatively, the correction to the a-enhanced pat- 
tern can be applied a posteriori: SP models com- 
puted with theoretical stars are used to differen- 
tially modify SP models based on the empirical li- 
brary. An application of this method was presented 



in ICervantes et al.l (|2007l ). using newly computed 
theoretical stars to correct SP models based on the 
MILES library (Vazdekis et al. in prep.). 

4. CONCLUDING REMARKS 

Considerable efforts have being applied in the 
modelling of SP models with variable chemical pat- 
terns. Following the pioneer efforts mentioned in 
§3.1, the models that combined fitting functions and 
sensitivity tables (§3.2) were a breakthrough, and 
they had an important impact on our knowledge 
of the abundances in galaxies derived by spectral 
indices. Recently, fully-theoretical models (§3.3) 
started to fill a gap which existed among the full 
spectral high-resolution SP models, providing ways 
of modelling a much larger number of observables 
and including spectral and evolutionary effects in a 
consistent way. And there are other methods be- 
ing developed nowadays, like the empirical correc- 
tions mentioned in §3.4. Each of the methods has 
its strong and weak points and we can expect fu- 
ture improvements in all of them. I think that there 
is not a 'best' approach. The user of stellar popula- 
tion models should choose the model family that best 
suits his/her application, always keeping in mind the 
weakness and strengths of each approach. But we 
can be certain of one thing: one way or another, we 
do need model stars in our galaxy models. 

Acknowledgements: P. Coelho thanks the organ- 
isers for their kind hospitality and financial support. 
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